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ABSTRACT 

Improved orbital parameters of the high-mass binary X-ray pulsar 4U 1538—52 have been derived from high 
count rate data obtained from the Rossi X-Ray Timing Explorer. Pulse-timing analysis yields an eccentricity of 
0.174 ± 0.015, a periastron at 64° ± 9°, and evidence of orbital decay with P orb IP„, „ = (-2.9 ± 2,1) x 10 

yr '• 

Subject headings: binaries: eclipsing — pulsars: individual (4U 1538-52) — X-rays: stars 


1. INTRODUCTION 

The accretion-powered high-mass X-ray pulsar 4U 1538—52 
(hereafter X1538) and its B0 I companion QV Nor comprise 
an eclipsing binary system exhibiting several interesting fea- 
tures that were revealed in previous observations by the Ginga 
X-ray observatory: the spectrum has a pulse-phase-dependent 
feature around 20 keV that was identified as an effect of the 
cyclotron resonance on the Compton scattering in the accretion 
columns of the neutron star (Clark et al. 1990); the X-ray image 
has a prominent grain-scattered halo that exhibits the effects 
of travel-time differences on its decay when the central source 
is eclipsed (Clark, Woo, & Nagase 1994); and pulse-timing 
analysis yielded evidence of orbital eccentricity with a value 
of 0.082 ± 0.047 (Corbet, Woo, & Nagase 1993). Optical ob- 
servations of the velocity curve of QV Nor, combined with 
previously determined orbital parameters of XI 538, yielded a 
mass of 19.8 ± 3.3 M @ for the companion star (Reynolds, 
Bell, & Hilditch 1992). 

This Letter describes new measurements of the orbital pa- 
rameters of XI 538 derived from pulse-timing analysis of high 
count rate data from the Rossi X-Ray Timing Explorer (RXTE). 
The results show a surprisingly large eccentricity and evidence 
of orbital decay. 

2. OBSERVATIONS AND DATA REDUCTION 

The observations occurred in 1997 during January 1 .93-5.69 
UT (MJD 50,449.93-50,453.69) with a total useful exposure 
of 139,723 s, extending from the binary orbital phase 0.93 over 
one complete orbital period. The RXTE observatory has been 
described by Bradt, Rothschild, & Swank (1993). It has five 
proportional counter arrays (PCAs) with a total sensitive area 
of 7000 cm 2 , a collimated field of view with a FWHM of 1°, 
and a nominal energy range extending from 2 to 60 keV. The 
raw data used in this analysis consist of the arrival time and 
pulse height of each event detected in the top two layers of 
the PCAs with a time resolution of 1 ms. The data were pro- 
cessed to obtain the count rates in 1 s intervals tagged with 
barycentrically corrected times. 

3. ORBIT ANALYSIS AND PULSE FREQUENCY 

Forty average pulse profiles were formed by epoch-folding 
the count rates into 80 bins of pulse phase. Each profile was 
derived from a nearly continuous segment of data spanning 
approximately 1% of the binary period. Pulse phases were 
computed for trial values of the intrinsic pulse frequency v, its 


rate of change v, and predicted delays in pulse arrival times 
derived by the solution of the Keplerian equation for trial values 
of the period P B , b , epoch T t/2 (time of mean longitude = 90°), 
projected semimajor axis a, sin i, eccentricity e, and periastron 
u of an elliptical orbit. A 10-term Fourier series was fitted to 
a template pulse derived from a continuous segment of data 
that was free of any substantial flare activity and extending 
from 120,000 to 134,000 s from the start of the observation. 
The Fourier representation of the template was cross-correlated 
with each of the 80 bin profiles, and the relative phase shift 
of each profile that was due to errors in the trial values was 
determined by a parabolic interpolation between the phase bin 
of maximum correlation and its immediate neighbors on either 
side. The predicted delay at the midtime of the segment was 
corrected according to these phase shifts in order to obtain a 
set of measured delays. A curve of predicted delay was then 
fitted by least squares to the set of measured delays in a two- 
step procedure by an adjustment of the frequency and orbital 
parameters with the Levenberg-Marquardt algorithm (Press et 
al. 1986). In both steps, P 0[h was fixed at the value determined 
as described below. In the first step, v was fixed at zero. In the 
second step, a x sin i was fixed at the fitted value from the first 
step, and i> was allowed to vary. This two-step procedure was 
necessary because a fit with both a t sin i and v free to vary 
was not adequately constrained by the data and yielded absurd 
values of v, v, and the orbital parameters. 

Four profiles with measured delays that deviated widely from 
the fitted curves of predicted delays were eliminated, two with 
low count rates that precluded accurate delay measurements 
from just before and just after the X-ray eclipse and two with 
large distortions in shape from intervals containing high flare 
intensities at orbital phases 0.17 and 0.42. 

Errors of the measured delays due to counting statistics were 
negligible compared with errors caused by variations in the 
intrinsic shapes of the profiles. To quantify these errors, a con- 
stant error o p was assigned to each of the delays and adjusted 
in magnitude to achieve = 1 in the second step of the fitting 
procedure. The adjusted value of o p is a measure of both the 
errors of the measured delays and the quality of the fit. 

A quadratic function of orbit number N was fitted to the set 
of all previous orbital epochs derived by pulse-timing analysis 
together with the present one. Mideclipse times estimated from 
observations of eclipse ingress and egress were not included 
because of their large uncertainties caused by the variable at- 
tenuation of the X-rays along the line of sight in circumsource 
matter near the times of eclipse transitions. The improved value 
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TABLE 1 

Parameters of the 4U 1538-52/QV Nor Binary System® 


Parameter Elliptical Orbit Circular Orbit 


^ sin i ( lt ' s > 56.6 ± 0.7 54.3 ± 0.6 

5,0450.206 ± 0.014 50,450.203 ± 0.009 

'W < s) 528.809 ± 0.008 528.754 ± 0.005 

s ^ (“0.4 ± 2.3) x 10 -13 (10.2 ± 1.7) x 10~ 13 

e 0.174 ± 0.015 

w (cleg) 64 ± 9 

X; (cr p = 2.27 s) J .00 3.29 

a 0 (MJD) b 47,221.488 ± 0.015 

fl > < da ys) b 3.728366 ± 0.000032 

a 2 (days)" (-5.5 ± 4.0) x 10“* 

b(yr“‘) (-2.9 ± 2.1) x I0” 6 


a Epoch = MJD 50,449.93400. 
b T wr2 { AO = a 0 + a , N 4 - a^N 2 . 


of the current orbital period computed from the fitted quadratic 
was used together with the other fitted parameters in a new set 

of epoch folds to sharpen the template and the average pulse 
profiles. 

This entire process was repeated several times until a final 
value of a p = 2.27 s was obtained. The resulting values of the 
orbital parameters are listed in Table 1 with 1 o errors derived 
from the diagonal elements of the error matrices. The final set 
of measured delays and the fitted curve are plotted against time 
from T sl&n = MJD 50,449.93400 in Figure 1 a. The deviations 
are plotted below. 

The orbital epochs are listed in Table 2, and their deviations 
from a fitted linear function of orbit number are displayed in 
Figure 2 together with the fitted quadratic curve. For the linear 
fit, Xy — 0.79. The quadratic fit has xl = 0.56 and gives 
PoJPor* = (-2.9 ± 2.1) X FT 6 yr" 1 . The low values of xl 
indicate that some of the errors attributed to the epochs were 
overestimated. 

Figure 1 b displays the delay curve and deviations for a cir- 



Fig. 1 — Observed delays in pulse arrival times and the fitted arrival time 
delay curves (a) for an elliptical orbit and ( b ) for a circular orbit. The fitted 
orbital parameters are listed in Table I. 


cular orbit (e = 0) fitted in the same manner to the same set 
of measured delays. With a p = 2.27 s, the circular fit had 

xl = 3.29. The fitted parameters for the circular orbit are listed 
in Table 1. 

4. DISCUSSION 

The systematic wave of the deviations of the measured delays 
from the circular fit and the superior quality of the elliptical 
fit are evidence of a substantial eccentricity of the XI 538 orbit. 
However, the value obtained is surprisingly large for a close 
binary system in which tidal friction is expected to cause rapid 
circularization. It is well outside the 95% confidence upper 
limit of e = 0.058 derived by Rubin et al. (1997) from ob- 
servations by the BATSE instrument on the Compton Gamma 
Ray Observatory and substantially larger than the value found 
by Corbet et al. (1993). The BATSE upper limit may be in 
error because of the difficulties inherent in the analysis of data 
with a very low signal-to-noise ratio recorded over many orbital 
periods during which fluctuations and the secular variation of 
the pulse frequency occurred. The BATSE analysis did yield 
a value of a x sin i = 53.5 ± 1.4 lt-s that is consistent with the 
value found here for the circular orbit fit. Nevertheless, the 
question arises as to whether the eccentricity found here as 
well, perhaps, as that found previously in the Ginga data could 
be artifacts of fluctuations in pulse frequency like those re- 
vealed in the BATSE study. 

In the BATSE observations, the pulse frequency was mea- 
sured at 16 day intervals over 1600 days. Rubin et al. (1997) 
found that the power density spectrum of the frequency deriv- 
ative fluctuations was approximately flat over the range of sam- 
pling frequencies from 5 x 10~ y to 3.6 x 10 -7 Hz s -1 with an 
rms noise level of S = (7.6 ± 1.6) x 10~ 21 (Hz s' 1 ) 2 Hz -1 . 
They found that the observed frequency history could be ac- 

TABLE 2 

Epochs Measured by Pulse-Timing Analysis 


Epoch 

Source Orbit Number (MJD) Error Reference 


OSOS -1128 43,015.800 0.100 Becker et al. 1977 

Tentna -457 45,517.660 0.050 Makishima et al. 1987 

Ginga 0 47,221.474 0.020 Corbet et al. 1993 

BATSE 370 48,600.979 0.027 Rubin et al. 1997 

478 49,003.629 0.022 Rubin et al. 1997 

584 49,398.855 0.029 Rubin et al. 1997 

691 49,797.781 0.022 Rubin et al. 1997 

RXTE 866 50,450.206 0.014 Present work 
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Fi G> 2— Differences between the epochs measured by pulse-timing analysis 
and the epochs predicted by a fitted linear function of orbit number. The curve 
is a plot of the differences between a quadratic function of orbit number fitted 
to the measured epochs and the linear prediction. The coefficients of the 
quadratic function are listed in Table 1 . 

counted for as a random walk in pulse frequency with a step 
size equal to a Gaussian variate with a zero mean and standard 
deviation ( SAt) y 2 , where At is the time interval between steps 
that they set equal to one orbital period, -1/4.3 of the shortest 
sampling interval of their frequency derivative measurements. 

We explored the possible effects on our determination of 
orbital parameters of a similar random walk in frequency, as- 
suming the same rms noise level to be valid for the even shorter 
time intervals between our delay measurements. Sets of 36 
artificial delays with random walk deviations were generated 
at the times of the observed delays. Specifically, to each exact 
delay computed at time t t for a given orbit, we added a deviation 

J — i 

AD, = X tty + e„ (1) 

j= l 

where 

&Dj = M ~ 0— impulse* / ^ 2, ... , 36. (2) 

Here represents the random measurement error computed as 
a Gaussian variate with zero mean and a standard deviation 
equal to 2.27 s, hD x = 0, hv x - 0, bv j>x is a Gaussian variate 
with zero mean and a standard deviation equal to [S(tj — 
tj-\)Y n , and P pulsc is the pulse period. 

One thousand sets of delays generated for a circular orbit 
with a x sin i = 56.6 lt-s were fitted with elliptical orbits as 
previously described. The rms frequency excursion was 
4.47 x 10“ 8 Hz, close to the expectation value of 4.44 x 
10“ 8 Hz for the time interval of 259,740 s between the first 
and last delay measurements. The rms delay excursion was 
3.60 s, close to the expectation value of 3.05 s if each frequency 


excursion had occurred at a constant rate of change. The mean 
values and rms deviations of the results were xt = 102 ± 
0.26, e = 0.033 ± 0.018, a x sin i = 56.65 ± 1.06 lt-s, 
v = (0.058 ± 0.744) x 1(T 13 Hz s“\ and T rl2 = MJD 
50,450.207 ± 0.015. The values of a) were distributed from 
0° to 360°, as expected. The largest e in the 1000 trials was 
0.10. The value of e derived from the real data is larger than 
its mean in the test sample by 7.8 times its rms deviation from 
the mean. Thus, it appears very unlikely that the eccentricity 
derived from the RXTE data is an artifact of random frequency 
derivative fluctuations, with the noise level found in the BATSE 
data. 

Previous measurements of change in the orbital period of 
X1538 have yielded values with 1 o error ranges that include 
zero. The value found here, (—2.9 ± 2.1) x 10 6 yr \ lies 
within the 95% confidence range of -3.9 < P ar JP OTh < 2.1 in 
units of 10 -6 yr 1 derived from the BATSE data (Rubin et al. 
1997) and is within a factor of 2 of the orbital decay rates that 
have been measured in Cen X-3 (Nagase et al. 1992), SMC 
X-l (Levine et al. 1993), and 4U 1700-37 (Rubin et al. 1996). 

The gross properties of X1538/QV Nor are similar to those 
of the binary pulsar SMC X-l/Sk 160 (see Levine et al. 1993). 
Both systems have a highly magnetized neutron star accreting 
matter from a massive B0 companion in an orbit with nearly 
the same period and rate of decay. However, X1538/QV Nor 
is clearly in an earlier stage of evolution. Whereas the orbit of 
SMC X-l is highly circular, with a 2 a upper limit of 4 x 
10~ 5 on its eccentricity, that of X1538 is substantially eccentric. 
And whereas SMC X-l draws its huge luminous power by 
accretion from a relatively stable precessing accretion disk 
(Wojdowski et al. 1998), X1538 evidently accretes primarily 
from the stellar wind of its companion and has a much lower 
X-ray luminosity. Driven by a steady inflow of angular mo- 
mentum from its accretion disk, SMC X-l has a high and 
monotonically increasing spin rate of 1 .4 Hz. Driven by inflows 
of angular momentum that flip in sign on timescales that may 
be as short as one orbital period or less, XI 538 spins up and 
down erratically around values close to 0.0019 Hz. The decay 
and circularization of the orbit and evolution of QV Nor will 
presumably bring a time when accretion occurs through a stable 
accretion disk fed by Roche lobe overflow, resulting in rapid 
monotonic spin-up like that observed in SMC X-l. 

The author thanks J. Woo for early collaborations in the 
planning and analysis of these observations and A. Levine, 
S. Rappaport, and D. Chakrabarty for valuable discussions. 
Thanks also to a referee for suggested improvements. Data for 
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the Goddard Space Flight Center. The help of the HEASARC 
staff in the data reduction is gratefully acknowledged. This 
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